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Abstract. The debate on the nature of the gamma-ray emission from young supernova 
remnants is still open. Ascribing such emission to hadronic rather than leptonic pro¬ 
cesses would provide an evidence for the acceleration of protons and nuclei, and this fact 
would fit with the very popular (but not proven) paradigm that supernova remnants are 
the sources of Galactic cosmic rays. Here, we discuss this issue with a particular focus 
on the best studied gamma-ray-bright supernova remnant: RX J1713.7-3946. 


1 Introduction 


The acceleration of protons at supernova remnant (SNR) shocks is accompanied by the emission of 
gamma-rays, resulting from the decay of neutral pions produced in interactions between the accel¬ 
erated cosmic rays (CR) and the ambient interstellar matter (see e.g. IT] and references therein). It 
is widely believed that the efficiency of energy conversion into CRs should be at the a; 10% level at 
SNR shocks, i.e. each SNR should convert Wcr ~ 10^** erg into CRs. The reason for that is phe¬ 
nomenological: it is the energy needed to sustain the observed intensity of CRs against their escape 
from the Galaxy 0. An estimate of the very high energy (VHE) gamma-ray emission resulting from 
the interactions between the CRs accelerated in a SNR and the ambien^as of density 4 x rigas (the 
factor of 4 accounts for the compression of the gas at the SNR shock) isQ 
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where t„o ss 1.6 x lO^/n^aj/cm^) yr is the energy loss time due to neutral pion production and d 
the distance to the SNR. In computing Eq. [^we implicitly assumed that the CR spectrum is cx 
and extends from = 1 GeV to E„,ax = 4 PeV, which is the energy of the CR knee. This is a 
large flux, well within the reach of current Cherenkov telescopes, unless the SNR is too far away 
or the ambient density is much lower than the average one in the Galaxy 0. A spectrum much 
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steeper than and/or a cutoff in the spectrum at an energy much lower than that of the knee would 
result is a dramatic suppression of the VHE gamma-ray flux, but these spectral feature would be in 
contradiction with the SNR paradigm for the origin of CRs. For this reason, the detection of SNRs in 
VHE gamma-rays was long awaited and considered a test for CR origin 13. However, the detection of 
several SNRs in VHE gamma-rays cannot be considered a proof of CR acceleration because leptonic 
processes (most notably inverse Compton scattering) could as well explain the observed emission H. 

The contribution to the gamma-ray flux from inverse Compton scattering of electrons in the cosmic 
microwave background radiation can be computed by recalling that photons of energy of a; 1 TeV are 
produced by electrons of energy ~ 20 TeV. At this energy the rate of production of TeV photons 
is Tj/! a! (5 X 10"^yr)“^ If the spectrum of electrons accelerated at a SNR shock is K^p times that 
of protons, then the contribution to the a; 1 TeV gamma-ray flux from inverse Compton scattering 
is ss KepWcR/2Tic\n(E^aJE„i„)(47Td^) ^ 2x eTg)(d/kpc)-^(Kep/10^'^) erg/cm^s. 

From which the ratio of leptonic to hadronic contribution to the VHE gamma-ray emission at 1 TeV 
can be deduced; 

fc,„(lTeV)»0.5(^)"(i) (2) 

which indicates that in normal interstellar medium conditions the leptonic contribution dominates un¬ 
less electrons are less than ^0.1% abundant than protons. In other words, inverse Compton scattering 
is much more efficient than proton-proton interactions in producing gamma rays and thus the leptonic 
emission from a SNR might well exceed the hadronic one even in the presence of a significant accel¬ 
eration of protons. Hence, the difficulty of ascribing unambiguously the VHE gamma-ray emission 
from SNRs to hadronic processes. 

The estimate given above is valid under the assumption that the shape of electron and proton spec¬ 
tra are identical. This is true only if the electron spectrum is not affected by energy loss processed 
In SNRs the main cooling channel for electrons is synchrotron radiation, which proceeds at a rate 
Tsyn ~ 6 X 10^(B/10 p.G)~^(Egl20 TeV)^' yr. This implies that in young SNRs of age of the 
order of few thousand years the multi-TeV spectrum of electrons is affected by synchrotron losses if 
the magnetic field is > 10 /uG. In this case, the ratio between leptonic and hadronic emission is sig¬ 
nificantly reduced with respect to the value computed in Eq. because a break forms in the electron 
spectrum at an energy £* that satisfies the condition Tage ~ Tsy„iE*). 

The spectrum of electrons accelerated at a SNR shock can thus be described as follows: at low 
energies, where both electrons and protons are unaffected by energy losses, the electron spectrum 
is identical in shape to the proton one, but is rescaled by the factor K^p which accounts for the dif¬ 
ferent acceleration efficiency of the two species. If the magnetic field strength is large enough, a 
break may appear in the electron spectrum at an energy E* ^ 1 (B/IQQ jjGy^iTage/kyry^ TeV. 
In gamma rays, the break in the inverse Compton spectrum would appear at a photon energy 
£* ss 4 (B/lOO iJ.Gy'^(Tage/kyr)~^ GeV. At energies larger than the break, the electron and gamma-ray 
spectra steepen by one power and half power of energy, respectively, and proceed until a maximum 
energy which can be computed by equating the energy loss time (determined by synchrotron radia¬ 
tion) to the acceleration time at the shock DB{Ee)lu^, where Db — {II3)Rlc is the Bohm diffusion 
coefficient, Rp the particle Larmor radius, c the speed of light, and Us the shock speed. This gives 
£max _ 20 (B/lOO /iG)“'^^(Ms/1000 km/s) TeV, which corresponds to a gamma-ray photon energy 
of a few TeV in the inverse Compton spectrum. Thus, a cutoff in the electron spectrum at energies 
< 30 TeV would prevent the production of VHE gamma rays through inverse Compton scattering, 
unless soft radiation fields other than the cosmic microwave radiation (e.g. infrared, optical ... ) are 
present. 

^Energy losses of protons have characteristic times well in excess than the SNR age and can be safely neglected. 
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Having identified all the main physical quantities that shape the hadronic and leptonic contribution 
to the gamma-ray spectrum of SNRs, we can proceed discussing a specific object, RX J1713.7-3946, 
in an attempt to understand the origin of its gamma-ray emission. 

2 The SNR RX J1713.7-3946 

RX J1713.7-3946 is the best studied SNR in VHE gamma-rays and the hadronic or leptonic nature 
of its emission is still debated (for a non-exhaustive list of references see ifS UTOll '). In the following 
we discuss this issue by listing the various arguments that have been brought forward in favor or a 
hadronic or leptonic origin of the emission. 

2.1 Arguments in favor of an hadronic gamma-ray emission 

2.1.1 Evidence for magnetic fieid ampiification 

According to shock acceleration theory, the acceleration of CRs at SNR shocks must be accompanied 
by the amplification of the magnetic field up to values of hundreds of milliGauss or even milliGauss. 
In order to have an effective amplification of the field, an efficient acceleration of CR hadrons is 
required IfTTlfT^ . This finding received strong observational support with the detection of thin X-ray 
synchrotron filaments from several SNR shocks. In an amplified field, electrons accelerated at shocks 
are cooled much more rapidly than the rate at which they are advected downstream of the shock. 
Thus, they radiate synchrotron X-rays in the close proximity of the shock and their emission appears 
as a very thin filament. The thinness of filaments allows to constrain the magnetic field strengths in 
the hundreds of microGauss range from a number of young SNRs (see ifTSll for a review). Though 
the presence of X-ray filaments in RX J1713.7-3946 was revealed by Chandra observations IfTTlI . the 
most spectacular evidence for magnetic field amplification in this object came from the observation, 
also made by Chandra, of fast (year-scale) variability of the synchrotron emission of some small (sub¬ 
parsec scale) X-ray knots El. If the variability is interpreted as the result of the fast cooling of freshly 
accelerated electrons, a magnetic field of a; 1 mG is inferred El- Alternatively, the rapid variability 
of the X-ray emission can be explained by a fluctuating magnetic field. This has been indeed proposed 
and discussed in ifTSlI . where the authors, in order to explain the variability, adopted a 100 pG field, 
which is lower than that inferred in M but still significantly amplified. 

It has to be stressed that the evidence for an amplified magnetic field is limited to a very small 
fraction of the SNR volume (filaments and knots) and thus it is difficult to reach strong conclusions 
on a global level. For example, it is unclear if CR protons are accelerated with large efficiencies in 
the whole shock surface or limited to a small fraction of it. Also, very little constrains are available 
for the strength of the magnetic field in the whole remnant. If the volume averaged value of the 
magnetic field is significantly larger =» 10 pG then a leptonic interpretation of the gamma-ray emission 
from RX J1713.7-3946 is problematic, because in this case electrons would radiate mainly through 
synchrotron radiation rather than inverse Compton scattering. On the other hand, if the magnetic 
field outside filaments and knots is at the level of ^ 10 pG then the leptonic interpretation would be 
definitely favored. 

To conclude, high resolution X-ray observations of SNRs, and in particular of RX J1713.7-3946, 
are generally interpreted as a very strong evidence for an efficient acceleration of CR hadrons from 
at least a part of SNR shocks. It is tempting to link such an efficient acceleration efficiency with a 
correspondingly high level of the hadronic gamma-rays from the remnants. However, the ambient gas 
density also plays a crucial role here (see Eq. [T]), and it can be easily shown that if the target density 
is as low as 0.1 cm“^ even an effective conversion of energy into CRs (=» 3 x 10^° erg in form of CR 
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protons) would not suffice to provide a significant hadronic contribution to the observed gamma-ray 
flux fl^ . 

2.1.2 Spatial correlation of gamma-ray emission and gas density 

The search for spatial correlations between gamma-ray brightness and gas column density is a pow¬ 
erful tool that can serve to discriminate amongst hadronic and leptonic emission. This is because a 
correlation with the ambient gas density is naturally expected in a hadronic scenario where CR pro¬ 
tons produce gamma rays in interactions with the ambient gas. Such a correlation is not granted (but 
at the same time not excluded!) in leptonic models. 

A spatial correlation of SNRs and dense gas is likely if the SNR progenitor is a supernova of type 
II, which is expected to explode inside a molecular cloud as a result of the very fast evolution of a 
very massive star El. Indeed, a type II progenitor for RX J1713.7-3946 seems to be suggested by 
observations and a correlation with dense molecular gas has been claimed E). The significance 
of the correlation has been found to increase as a result of detailed studies that considered the total gas 
density (atomic plus molecular) in the region ll20l . Also, the ambient gas within the SNR shell appears 
to be structured in small (sub-parsec) and dense (up to at least 10"^ cm“^) clumps El- The presence 
of a significant amount of gas in the SNR shell might of course provide a target for CR proton-proton 
interactions and enhance the hadronic contribution to the gamma-ray emission GO). 

2.2 Arguments in favor of a leptonic gamma-ray emission 

2.2.1 Spectral shape of the gamma-ray emission 

The spectral shape of the gamma-ray emission from RX J1713.7-3946 is probably the most popular 
argument presented in favor of a leptonic origin of the emission. As described in the Introduction, 
if the magnetic field is amplified, two features are expected in the leptonic gamma-ray spectrum of 
SNRs; i) a break at energy £* due to the fact that particles above a given energy cool via synchrotron 
emission on a time scale shorter than the age of the remnant, and ii) a cutoff at an energy which 
represents the maximum energy a particle can be accelerated to before losing energy via synchrotron 
radiation. If we consider an injection spectrum of electrons proportional to then the electron 
(gamma-ray) spectrum is oc E^^ (oc below the break and oc E^^ (oc E^"^) above it. 

The picture described above changes if the magnetic field at the SNR is weak. This is because 
both the break and the cutoff energies increase for weaker magnetic fields, but in two different ways, 
i.e., £* oc and oc Z? implying that a magnetic field strength exists below which the break 
in the spectrum disappears. Such a field is ^ 10 pG (see Introduction), and for field strength smaller 
or equal to this value the electron (gamma-ray) spectrum becomes a simple power law in energy with 
index -2 (-1.5). 

Indeed, a weak magnetic field is required for the leptonic interpretation to be viable. Synchrotron 
X-ray emission has been detected from RX J1713.7-3946 at a flux level of E^Fx ~ 10“'® erg/cm^s, 
while the VHE gamma-ray emission E^Fy lays roughly one order of magnitude below that. If both 
X-rays and gamma rays are produced by the same electrons through synchrotron and inverse Compton 
scattering in the cosmic microwave background, respectively, then their flux ratio should reflect the 
ratio between the energy density in magnetic field ub - B^f^n ^ 2.5 (B/10 pG)^ eV/cm^ and in the 
cosmic microwave background ojcmb ~ 0.25 eV/cm^: 
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which implies that a field strength of a; 10 jiG is required to fit data. 

To summarize, a simple, one-zone model leptonic interpretation of the spectrum of RX J1713.7- 
3946 requires a weak magnetic field, and a rising spectral energy distribution of gamma-ray photons 
below the cutoff: E^Fy oc Ey^. On the other hand, if the magnetic field is » 10 yuG, an hadronic 
interpretation seems favored, and for a typical oc injection spectrum of protons one would expect 
a flat spectral energy distribution of gamma-rays from neutral pion decay: EyFy oc Ey. This is the 
reason why the detection by Fermi of a rising spectral energy distribution in the GeV domain from 
RX J1713.7-3946 (see Ifl^ and Fig.[T]) has been interpreted as a strong evidence for the fact that the 
gamma-ray emission is indeed leptonic. Following the same reasoning, the much steeper gamma-ray 
spectrum detected from the Tycho SNR has been interpreted as an evidence for hadronic emission 

ESI US. 

In fact, it has been soon realized that a one-zone leptonic model fails to obtain a good fit to 
the gamma-ray data for RX J1713.7-3946. The reason is that the expected spectrum has a peak much 
narrower than the observed one (see e.g. Fig. 3 in |j4]|). This problem can be solved either by assuming 
that two populations of relativistic electrons exist (e.g. accelerated at the forward and reverse shock, 
respectively?) ll9]|24l, or by invoking the presence of a strong infrared radiation field in the SNR 
that would provide a target for inverse Compton interactions distinct than the cosmic microwave 
background Ea. However, such a strong additional radiation field seems quite difficult to be justified 
(see e.g. 0). 

2.2.2 Absence of thermal X-ray emission 

The absence of thermal Bremsstrahlung continuum emission in the X-ray spectrum of RX J1713.7- 
3946 was claimed to be an evidence against the hadronic origin of the gamma rays ll26l . In order to 
have effective CR hadronic interactions the ambient gas density must not be too low, since the hadronic 
emission scales as the target gas density. On the other hand, if this is the case one should also expect 
a powerful X-ray thermal emission from the gas heated by the SNR shock (thermal emission scales 
as the gas density squared). In fact, it has been shown that the problem comes mainly from the non 
detection of X-ray thermal lines, rather than continuum emission. Lines should overcome the observed 
synchrotron emission of RX J1713.7-3946 unless the ambient density is very low, < 0.05 cm“^. Since 
no thermal emission is observed (neither continuum nor lines) hadronic models are challenged ||8]. 

The hadronic interpretation would still be viable if the gas is not heated to the high temperatures 
required for the production of X-rays. This could happen if the shock converts a large fraction (signif¬ 
icantly above « 10%) of the incoming kinetic energy into CRs, since in this case not enough energy 
would be left to heat the gas (e.g. EIHlol). However, recent simulations of the injection phase of 
particles at shocks suggest for the CR acceleration efficiency a value of a; 10% ED, implying that 
significan gas heating must occur. Another interesting possibility to avoid gas heating involves the 
presence of dense clumps in the medium surrounding the SNR and will be discussed below. 

2.3 A hadronic model for RX J1713.7-3946 

It has been argued in Q and E 2 I that the origin of the gamma-ray emission from RX J1713.7-3946 
may well be hadronic if the SNR expands in a clumpy medium. The simulations performed in 0^ 
suggests that this is a quite natural situation if the SNR progenitor is a type II supernova exploding in 
a molecular cloud, as RX J1713.7-3946 is believed to be IfTBl . The powerful wind from the progenitor 
star would sweep away the dense gas creating a cavity characterized by a gas density of the order of 
ss 10“^ cm"^. However, the densest cores (of density above =» 10^ cm“^) in the molecular cloud would 
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Figure 1. Hadronic model for RX J1713.7-3946, under the 
assumption that the remnant originated from a type II 
supernova explosion in a molecular cloud Oni. The SNR is 
assumed to expand in a rarefied (~ 2 x 10“^ cm“^) cavity 
inflated by the wind of the progenitor star. The wind swept 
away all the matter but the densest parts (> 10^ cm“^), 
concentrated in small (~ 0.1 pc) clumps. Such dense clumps 
survive intact the passage of the SNR shock. CRs accelerated 
at the SNR shock diffusively penetrate into the clumps. Bohm 
diffusion in a x: 100 /rG field is assumed to operate at the 
interface between the clumps and the diffuse medium in the 
SNR. The solid line represents the predicted emission from 
the clumps, while the dashed line that from the diffuse gas 
inside the SNR. Here, an explosion energy of 10^' ergs, a CR 
acceleration efficiency of 10% and a spectrum of accelerated 
CRs at the shock oc ^ have been assumed. To match the 
data, the total mass in the clumps has to be x: 550 Mq. Data 
points are from FERMI (TD and HESS (Ml- 


survive into the cavity. Remarkably, such clumps would also survive the SNR shock passage. This 
happens because a SNR shock of velocity Us that hits a dense clump generates a shock into the clump 
itself, with a velocity Uc ~ uja, where a ^ (10^cm“^/10“^cm“^)'^^ a; 300 is the square root of the 
ratio between the clump to inter clump density (33]. Thus, the SNR shock is virtually stalled at the 
clump border and the clump remains unshocked. 

Remarkably, this fact solves one of the problems of hadronic models. The absence of thermal 
emission can be explained if the gas mass within the SNR is dominated by the contribution of dense 
clumps. Given that clumps remain unshocked and thus cold no appreciable X-ray thermal emission is 
expected. 

The presence of clumps can also solve the problem related to the observed gamma-ray spectrum 
of RX J1713.7-3946, which is generally considered too hard to be explained by hadronic interactions. 
According to simulations, the interaction between the SNR and the clumps would lead to a strong 
amplification of the turbulent magnetic field at the interface between the clumps and the diffuse gas 
(321 . Such a turbulent magnetic field would prevent low energy CRs (accelerated at the SNR shock 
and confined within the SNR shell) to penetrate the clumps over a time scale shorter than the SNR age 
(=s 1600 yr). On the other hand, high energy CRs would penetrate the clumps, their transport being 
characterized by a faster diffusion coefficient. The outcome of this is that the CR spectrum inside 
clumps is expected to be much harder than that of CRs in the SNR shell. If the total mass of the gas 
in the shell is dominated by the clumps, then the gamma-ray spectrum observed from the SNR would 
be hard, as the underlying CR spectrum inside clumps. 

Following these arguments, a detailed model for RX J1713.7-3946 has been developed in m, 
demonstrating that an excellent fit to the data can be obtained with a hadronic model (see Fig. [1]). 


3 Conclusions 

The origin of the gamma-ray emission from RX J1713.7-3946 is not well established yet. A leptonic 
scenario requires a low magnetic field strength (=» 10 fiG), and either two populations of relativistic 
electrons or two populations of target soft photons for inverse Compton interactions. On the other 
hand, a hadronic scenario implies a much stronger magnetic field and requires the presence of a 
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Structured ambient medium with dense clumps embedded in a tenuous diffuse gas. The detection or 
non detection of neutrinos will allow to discriminate amongst these two scenarios. For the hadronic 
scenario the detection of neutrino is well within the reach of a km^-scale detector located in the 
northern hemisphere 1351 . 

Also the detection of radiation from beyond the SNR shell would provide strong support to the 
hadronic scenario. Such a radiation is expected from the vicinity of the shell as the result of the 
interactions of escaping CR protons in the ambient gas IHEI. For RX J1713.7-3946 this would 
result in the presence of a faint diffuse gamma-ray emission surrounding the SNR, of brightness at the 
limits of detectability of current Cherenkov instruments, but within the reach of the future Cherenkov 
Telescope Array ll^ . Remarkably, also X-ray telescopes like NuSTAR and Astro-H can contribute to 
the search of escaping CR protons by detecting the synchrotron emission produced by the secondary 
electrons produced by runaway CRs in interactions with the ambient gas ED. 
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